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Abstract 

We calculate the full one-loop electroweak radiative corrections, of 0{a'^as), to 
the cross section of single Z-boson inclusive hadroproduction at finite transverse 
momentum (pt)- This includes the 0(a) corrections to Z + j production, the 
0{as} corrections to Z -|- 7 production, and certain QCD-electroweak interference 
contributions involving a single quark trace. We recover the QCD and purely weak 
corrections and study the QED corrections and the QCD-electroweak interference 
contributions for the first time. We also consider direct and resolved photoproduc¬ 
tion in elastic and inelastic scattering. We present pT and rapidity distributions 
for the experimental conditions at the Fermilab Tevatron and the CERN LHC and 
assess the significance of the various contributions. 

PACS numbers: 12.15.Lk, 12.38.Bx, 13.85.Qk, 14.70.Hp 
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1 Introduction 


The study of single electroweak-gauge-boson hadroproduction, via the so-called Drell- 
Yan process, has a long history, starting from the discovery of the W [1] and Z [2] 
bosons at the CERN Super Proton Synchrotron (SPS) more than three decades ago, 
which marked a breakthrough for the Standard Model (SM). These processes remain 
to be of paramount importance also at modern hadron colliders, such as the Fermilab 
Tevatron and the CERN LHC. On the one hand, they have large cross sections and clean 
decay signatures in the detectors. This renders them particularly useful for calibrating 
and monitoring the luminosities at hadron colliders, which affects all other cross section 
measurements performed there as well. By the same token, they are sensitive probes of 
the parton density functions (PDFs), in partiular of those of the quarks and antiquarks. 
On the other hand, singly produced W and Z bosons form important backgrounds for 
searches of new physics beyond the SM, such as anomalous couplings, extra vector bosons, 
etc. 

To achieve an adequate theoretical description, radiative corrections, both of QCD 
and electroweak type, must be taken into account. As for the total cross sections of single 
W- and Z-boson hadroproduction, the next-to-leading-order (NLO) [3] and next-to-next- 
to-leading-order (NNLO) [1] QCD corrections were calculated a long time ago, and also 
partial results at next-to-next-to-next-to-leading-order are available [S] . These corrections 
are of relative orders 0(ag) with n = 1, 2, 3, respectively, in the strong-coupling constant 
ag. The one-loop electroweak corrections, of relative order 0(a) in Sommerfeld’s fine- 
structure constant a, were studied for the W boson in Ref. [6] and for the Z boson in 
Ref. [7]. At the mixed two-loop order O(aag), the corrections to the qqZ form factor were 
evaluated for light quarks q b,t in Ref. [8] using the techniques developed in Ref. [9], 
and the treatment of the nonfactorizable corrections in the resonance region was discussed 
in Ref. [TO] . 

In order for the W and Z bosons to acquire finite transverse momenta (pr), they 
must be produced in association with additional particles or hadron jets (j). The QCD 
corrections to the pt distributions of W- and Z-boson inclusive hadroproduction were 
computed at NLO uma and partly at NNLO [13]. The 0{a) corrections were investi¬ 
gated for the W boson in Refs. gins] and for the Z boson in Refs. nan] - Specihcally, 
in Ref. [H], the electroweak 0{a) corrections to the 0{aas) partonic subprocesses of 
hP-boson production were calculated imposing a minimum-transverse-momentum cut on 
outgoing gluons to prevent soft-gluon singularities. However, this cut was not applied to 
outgoing quarks and antiquarks as well, which renders it impractical at the hadron level, 
where gluon and light-quark jets are hard to distinguish on an event-by-event basis. Such 
a cut is also problematic from the conceptual point of view because, as a matter of prin¬ 
ciple, a collinear gluon-photon system cannot be distinguished from a single gluon with 
the same momentum. In Ref. [151 , the results of Ref. [2] were conhrmed, but the soft- 
gluon singularities were properly eliminated by including also the 0{as) corrections to 
the 0{a‘^) partonic subprocesses. Furthermore, the 0{a^) contributions due to direct and 
resolved photoproduction by elastic and inelastic scattering off the incoming (anti)proton 
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were taken into account in Ref. |15j . 

In contrast to the charged-current case, the separation of the electroweak 0{a) correc¬ 
tions to the neutral-current Drell-Yan process into an electromagnetic and a weak part 
is meaningful with regard to infrared (IR) and ultraviolet (UV) hniteness and gauge in¬ 
dependence. In Refs. [inin], the purely weak 0{a) corrections to the 0{aas) partonic 
subprocess qq —)■ Zg and its crossed versions were computed. They make up an important 
subset of the contributions of absolute order 0{a‘^as) to the inclusive hadroproduction 
of hnite-pT Z bosons. It is the purpose of this work to complete our knowledge of these 
contributions, which have several sources, and to check the results presented in Ref. HZI. 
To start with, we need to complement the purely weak 0(a) corrections to the 0{aas) 
partonic subprocess qq —?■ Zg and its crossed versions by the QED ones, which have vir¬ 
tual and real parts. The 0{a^as) partonic subprocesses that we are then led to consider 
include qq —?■ Zgj. As in the charged-current case [15], we thus inevitably encounter a 
soft-gluon singularity. To cancel it, we need to also include the 0{as) QCD corrections to 
the C>(a^) partonic subprocess qq —?■ Z'j. Furthermore, 0{a^as) contributions may also 
arise from interferences of 0{a^^'^as) and Feynman diagrams yielding a single 

Dirac spinor trace with nonvanishing color factor trT“T“ = NcCp. This happens for 
the partonic subprocess qq ^ Zqq when a diagram involving a virtual gluon in the s {t) 
channel is connected with a diagram involving a photon or a Z boson in the t (s) channel, 
and for the subprocesses qq ^ Zqq and qq ^ Zqq when a gluon-exchange diagram and 
a photon/Z-boson-exchange diagram are connected with twisted quark lines. For com¬ 
pleteness, we also recalculate the 0(as) QCD corrections to the inclusive hadroproduction 
of hnite-pr Z bosons [TTIIT^ and thus recover the analytic results specihed in Ref. [ 12 ] 
apart from a few misprints that we correct. Finally, we also study the leading-order (LO) 
photon-induced subprocesses, of order O(a^), which, after convolution with the photon 
PDFs of the (anti)proton, yield contributions of absolute order 0{a^). As in Refs. [T^IITT] . 
we list full analytic results in a compact form. 

Our goal is to study the inclusive hadroproduction of single Z bosons with hnite 
values of pt- For the experimental analysis, this implies that all events with at least one 
identihed Z boson are selected and sampled in bins of one or more kinematic variables 
exclusively pertaining to the Z boson, such as pt and rapidity y. If there is more than 
one identihed Z boson in such an event, then each of them generates one entry in the 
considered histogram. There is no need to identify particles of other species or jets that are 
produced in association with the Z bosons. If such additional experimental information 
is available, it is nevertheless ignored. Samples of events with at least one identihed Z 
boson may, of course, also be analyzed more exclusively. For instance, one may study 
the production of a large-pT Z boson in association with a jet or a prompt photon. The 
separation of Z -|- j and Z -|- 7 events is efficiently achieved by means of the procedures 
elaborated in studies of the photon fragmentation function [18] and of photon isolation 
|19] . The contributions from Z + X hnal states, in which the system X contains a heavy 
particle, e.g. a W, Z, or Higgs boson or a top quark, are greatly suppressed and not 
considered here. 

The 0(a) corrections to + -f-j, and w + j inclusive hadroproduction were 
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considered in Refs. [201I21II22] . respectively. The theoretical stndy of such final states is 
closer to the experimental situation, as it does not rely on the identihcation of the W and 
Z bosons and the reconstruction of their four-momenta. The latter two procedures have 
been routinely applied in experimental data analyses ever since the discovery of the W 
and Z bosons at the SPS in 1983. They are, of course, subject to certain experimental 
errors, which are, however, quite small for the gold-plated Z e~^e~ and Z 
decay modes of relevance here. By the same token, the numerical results presented in 
Ref. [21] do not allow one to extract the px distribution of the Z boson and thus cannot 
be usefully compared with the results obtained in Refs. [Mnz] and here. This would 
require kinematic cuts to reduce the contributions due to the nonresonant parts of the 
scattering amplitudes in Ref. [21] in analogy to the experimental acceptance cuts, e.g. 
the one confining the invariant mass Mu of the l^l~ pair to an appropriately narrow 
interval about Mz- Apart from that, in Ref. [21], the cross sections were not presented 
as distributions in the transverse momentum of the l^l~ pair, which could be identified 
with the kinematic variable px of the Z boson for the sake of a comparison with the 
results obtained in Refs. [MHZ] and here. On the other hand, the Mu distributions 
shown in Figs. 7 and 8 of Ref. [21] exhibit a rapid fall-off at the shoulders of the peak 
at Mu = Mz indicating that the narrow-width approximation adopted in Refs. [MUni 
and here is quite appropriate for the Tevatron and the LHC. The 0{a'^as) interference 
contributions mentioned above were neglected in Refs. [2IE2] appealing to the observation 
that similar contributions were found to be numerically small in Ref. izg. Recently, 0{a) 
corrections were also calculated for the hadropro duct ion of the hnal states + 2j [23] 
and -|- nj with n = 1, 2, 3 [2^ . 

This paper is organized as follows. In Sec. [21 we explain our analytic calculations. In 
Sec. El we present our numerical results. Our conclusions are contained in Sec. jH Our 
analytic results are listed in Appendices OO 


2 Analytic results 


We consider the inclusive production of a Z boson in the collision of two hadrons hi and 

h2, 

hi{Pi) + h 2 {P 2 ) ^ Z{q) + X, (1) 


where the four-momenta are indicated in parentheses and X collectively denotes the 
residual particles in the final state. We take the Z boson to be on mass shell, = M|, 
neglect the hadron masses, Pf = P 2 = 0, and define the hadronic Mandelstam variables 
as 

S = iPi + P 2 f, T={Pi-qf, U={P 2 -q)\ Q^ = q^. (2) 

In the center-of-mass frame, we write q^ = (g°,qT) 7'^)5 where qy is the transverse mo¬ 
mentum, and dehne qx = [qnl and the rapidity y = (l/2)ln[(g° -|- — q^)]. Using 

Eq. ([21), we have 


Qx 


TU-Q^{S + T + U-Q^) 
S 



U-Q^ 


(3) 
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We work in the collinear parton model of QCD with nj = 5 massless quark flavors. 
We write the partonic subprocesses that contribute to the hadronic reaction in Eq. ([1]) 
generically as 

iiPi) + jiP2) ^ Z{q) + X{px), (4) 

where pi = XiPi with i = 1,2. The partonic Mandelstam variables are dehned as 

s = {pi+ P 2 )^ t = {pi- qf, u = {p 2 - qf, S 2 = Px, (5) 

and satisfy 

s + t M + S 2 . (6) 

The hadronic and partonic Mandelstam variables are related as follows: 

s = xiX 2 S, t = xi{T - Q‘^) + Q'^, u = X 2 {U - Q‘^) + , 

S2 = X 1 X 2 S + xi(T — Q^) + X2(U — . (7) 

The differential cross section for reaction ([T]) may be evaluated according to 

= Y1 j ^^i(^^2fi/hA^i,pl')fj/h2{x2,pl)^^^^{xiPi,X2P2,pl), ( 8 ) 

x ij 

where the sum runs over all the partons i and j; fi/hix,pp) is the PDF of parton i in 
hadron h] and fip is the factorization scale. The partonic cross sections dcTij / {dt du) may 
be computed in perturbation theory as double series in ag and a. Apart from the Feynman 
rules of QCD, we also need those for the couplings of a quark q to the photon 7 and the 
Z boson. They are given by the vertices ieQq'f^ and ze 7 ^(ng — j^aq), respectively, where 
Qq is the electric charge of q in units of the positron charge e = x/dvra, 

^ ^ Jg - 2Qq sin^ Oy, ^ ^ I 3 

2 sin 6 *^ cos 6 ^ ' '^2 sin 6 *^ cos 6 ^ 

are its vector and axial vector couplings to the Z boson, I 3 is its third component of weak 
isospin, and 6 ^ is the weak mixing angle. 

In the remainder of this section, we list the relevant partonic subprocesses with the 
contributing Feynman diagrams and outline the computation and its organization. At 
LO, we consider the 2 —)■ 2 subprocesses 


q + q^Z + g, (10) 

g + g^Z + 7, (11) 

q + g^Z + q, (12) 

g + 7^Z + g, (13) 


where it is understood that q may also be an antiquark, in which case g is a quark. They 
are mediated by the Feynman diagrams depicted in Fig. [H At LO, subprocesses flTOll and 
flT^ are of O(aas), and suprocesses flTTll and flT^ are of 0 {a‘^). 
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Figure 1; Tree-level diagrams contributing to the partonic subprocesses q + q ^ Z + g 
and q + g^Z + q at 0{aas) and to the partonic subprocesses q + q —)■ Z + 'j and 
q + 'j^Z + qat 


We include the 0{as) QCD corrections to subprocesses (|T0l) - (lT^ and the 0{a) elec- 
troweak corrections to subprocesses ca and (IT^ . while we treat subprocess (IT^ only 
at LO because of the additional 0{a) suppression due to the photon emission by the 
incoming hadrons. The virtual QCD corrections to subprocesses oa and oa and the 
virtual QED corrections to subprocess flTOl) arise from the Feynman diagrams in Fig. [2l 
The virtual QCD and QED corrections to subprocess flT^ are obtained by appropriately 
crossing external legs. The virtual weak corrections to subprocess ca arise from the 
Feynman diagrams in Fig. [3l and those to subprocess ffT2|) emerge by crossing. 

The real QCD and QED corrections arise from the 2 —)■ 3 subprocesses 


q + q^Z + q + q, 

( 14 ) 

q + q^Z + q+q', 

( 15 ) 

q + q^ Z + g + g, 

( 16 ) 

q + q^ Z + g + 'y, 

( 17 ) 

q + q^Z + q + q, 

( 18 ) 

q + q'^Z + q + q', 

( 19 ) 

q + g^ Z + q + g, 

( 20 ) 

q + g^ Z + q + 'y, 

( 21 ) 

g + g^Z + q + q, 

( 22 ) 


where q' ^ q. Subprocesses ffTT|) and flTH]) contribute both at and 0{a'^as), 

subprocesses flTSD . flTBD . flT^ . (120]) . and contribute at 0{aa‘^), and subprocesses flTTD 
and (l2T]) contribute at 0{a‘^as). The tree-level diagrams contributing to subprocess ffldl) 
are shown in Fig. 01 those contributing to subprocesses flT^ and flTTl) in Fig. [5l and those 
contributing to subprocesses ffTSj) and flT^ in Fig. [6l 

As already mentioned in Sec. [H the 0{a‘^as) contributions to subprocesses flT^ and 
flT^ are generated by interferences of 2 —)■ 3 tree-level diagrams with a virtual gluon in one 
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Figure 2: QCD and QED one-loop diagrams contributing to the partonic subprocess 
q + q^Z + g at O(aa^) (i?i = B 2 = g) and 0 {a‘^as) (diagrams (a l)-(a 10 ) with i?i = 7 
and B 2 = g) and to the partonic subprocess q + q^Z + '-^ at 0(a‘^as) (diagrams (a 1 )- 
(a 10) with Bi = g and B 2 = j). The QCD and QED one-loop diagrams contributing 
to the partonic subprocess q + g —)■ Z + q at 0{aa'^) and 0(a‘^as) are obtained by 
appropriately crossing external legs. 


factor and a virtual photon or Z boson in the other one in such a way that one closed quark 
line is formed yielding a nonvanishing color factor, tr T“T“ = NcCp, as indicated in Fig. [71 
In the case of subprocess (ITTll . this is achieved when the gluon is in the s channel and the 
photon or Z boson is in the t channel or vice versa. In the case of subprocess flT 8 |) . this is 
achieved by twisting the quark lines in one of the interfering diagrams. These interference 
contributions exhaust the 0{a'^as) corrections to subprocesses (IT4)) and (ITSU and are thus 
hnite and gauge-independent by themselves. On the other hand, interferences leading to 
two closed quark lines are nullified by (trT“)^ = 0. This explains, why subprocesses ffTHll 
and flT^ do not receive 0{a‘^as) contributions. Obviously, these types of mixed QCD- 
QED corrections may not be obtained by merely manipulating coupling constants and 
color factors as is often the case for pure QED corrections. 

We compute the full 0{aal) and 0{a‘^as) corrections to the cross section of the 
hadronic process ([I]) according to Eq. ([H]) by including all the partonic subprocesses (ITOll - 
(122|) . We regularize both the UV and IR divergences using dimensional regularization 
in D = 4 — 2e space-time dimensions. The UV divergences arise from the 2 —)■ 2 one- 
loop diagrams and are removed by renormalizing the coupling constants, masses, and 
wave functions in the respective 2 —)■ 2 tree-level diagrams. We renormalize as and a 
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according to the modified minimal-subtraction (MS) scheme and employ the electroweak 
on-shell renormalization scheme otherwise. In particular, we define 6 ^ in terms of the 
pole masses as cos^^ = Mw/Mz- The IR divergences, both of soft and collinear types, 
are generated by 2 ^ 2 one-loop and 2^3 tree-level diagrams. The soft and collinear 
divergences related to the final states are canceled by integrating over the kinematic 
degrees of freedom of the three-particle phase space that are related to the systems X 
and combining the outcome with the virtual corrections. Specifically, the three-particle 
kinematics turns into the two-particle one by taking the limit S 2 —)■ 0. This is implemented 
in dimensional regularization using the relationship 

^ (-M .r - f .r+■■■) + (i) ^ 


where is the maximum value of S 2 allowed for given values of pt and y and the plus 
distributions are defined for smooth test functions /(S 2 ) as 




ds2 — [f{,S2) -/(O)], 
S 2 


-glllclX - 

Jo ^2 


( 24 ) 


There remain collinear divergences related to the initial states, which are universal 
and are absorbed into the bare PDFs so to render them finite. In the MS factorization 
scheme, this PDF renormalization is implemented in the QCD sector as 



6 ij6{l - y) 



^ F(1 - e) 
2n F(1 - 2e) 


Pijiy) + ■ ■ ■ 


(25) 


where Pij{y) are the j ^ i splitting functions. In the one-loop approximation of QCD, 
the latter read [25] 


Pqq{y) = Cf 


+ 2 I -1-9 

1 + (1 - yf 


y 

PmM = (yCx - jrn, ) 5(1 -y) + 2 Ca 

PM = riy + (1 - 9)^1, 


1 - 2 / 


+ - - 2 + y{l - y) 

y 


(26) 


where Cp = (A(? — l)/{2Nc) = 4/3, Ca = Nc = 3, T = 1/2, and n/ = 5 is the number 
of active quark flavors. For simplicity, we adopt the MS factorization scheme also for 
the QED sector. The appropriate counterparts of Eqs. fl25l) and fl26l) are obtained by 
substituting as M Cp —^ 1, and Ca —^ 0. Of course, the PDFs to be used in the 

















numerical analysis must then be implemented with the MS factorization scheme as well, 
both in their QCD and QED sectors. While the MS factorization scheme is now common 
standard for the QCD sector, alternative choices have been advocated for the QED sector, 
e.g. in connection with neutrino-nucleus deep-inelastic scattering (DIS), a DIS-like choice 


To exploit Eq. ([23]), it is useful to introduce S 2 as an integration variable in Eq. (| 8 |), 
in lieu of X 2 , say. This leads to 


da 


dq^ dy 


where 


mm _ 
0/1 — 


1 '=>2 

^ ./ / 


^min 

Xi 


-u 


dso 


xiS + U-Q- 


' fi/hi (^1) fj/h2 (^2) hi^) 


^ iXiPi,X2P2,yl), 


dt du 


(27) 


^ + T-Q 2 ’ 


s^-- = U + xi{S + T-Q^), X2 = 


S2 - - xi{T - Q^) 


xiS + U-Q^ 


T = Q^-e-y^S{Q'^ + ql), U = ^S{Q'^ + q^). 


(28) 


At this point, we compare our analytic results with the literature [T211T7] . The NLO 
QCD corrections, of relative order 0{as), due to the virtual contributions from subpro¬ 
cesses (na and inai, and the real contributions from subprocesses (|T4|) - (lT6l) . (|T 8 j) - (j 20 j) , 
and (| 2 ^ are listed in Ref. [ 12 ]. Apart from some misprintsji] we hnd agreement with 
Ref. [I 2 ]. The weak 0(a) corrections to subprocess (ITOD are listed in Ref. [T7]. In 
Ref. [I^, collinear divergences arising from box diagrams in intermediate steps of the 
calculation are regularized by introducing an inhnitesimal quark mass A, while we employ 
dimensional regularization. The A-dependent one-loop scalar box integrals J 12 , T 13 , and 
Ji 4 in Eq. (39) of Ref. [1^ may be conveniently converted to dimensional regularization 
using the results of Ref. [27]. In Ref. ca. the renormalization is performed both in the 
MS scheme and in the on-shell scheme implemented with some running hne-structure 
constant as explained in Eqs. (49) and (50) of Ref. [IT], which differs from the pure 
MS dehnition. Specihcally, in the counterterm of the electromagnetic coupling constant, 
the photon self-energy, which appears there with argument = g in the pure on-shell 
scheme, is split into the fermionic and bosonic parts, and the argument of the former is 
shifted to q^ = M|. While the latter construction is well dehned at one loop, it becomes 
ambiguous at higher orders because of the required separation of fermionic and bosonic 
contributions, and it is bound to render the running hne-structure constant thus dehned 
gauge dependent. By contrast, we work in a hybrid renormalization scheme, which uses 
the pure MS dehnition of a, but the electroweak on-shell scheme otherwise. Taking these 
conceptual diherences into account, we fully agree with Ref. HZI. 


^In Eq. (2.12) of Ref. [H], i?|'^(s, t, m, should be replaced with 
and C'|'^(s, t, u, Q^) with t, u, Q'^). 


t, u, Q^) + C^^{s, t, u, Q"^) 
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In this paper, we only list those analytic results that may not be found in the previous 
literature. Specihcally, we consider 0{a'^as) contributions to q+q Z+X, q+g —)■ Z+X, 
and q + q —)■ Z + X in Appendices lAl iBl andICl respectively. In the case oiq + q —)■ Z + X, 
this includes the virtual QED corrections to subprocess (BO, the virtual QCD corrections 
to subprocess ca. the real corrections from subprocess on, and the above-mentioned 
interference contributions from subprocess flT^ . In the case of q+g —)■ Z+X, this includes 
the virtual QED corrections to subprocess flT^ and the real corrections from subprocess 
fl2T]) . In the case oi q + q Z + X, this includes the above-mentioned interference 
contributions from subprocess (lT8l) . 

As already mentioned in Sec. [H we also include the LO contributions from photo¬ 
production. Incoming photons can participate in the hard scattering either directly or 
indirectly, i.e. through their quark and gluon content. The contributions from direct and 
resolved photoproduction are formally of the same order in the perturbative expansion. 
This may be understood by observing that the PDFs of the photon have a leading be¬ 
havior proportional to aln(p|./AQ(^j 3 ) oc a/as{+%), where Aqcd is the asymptotic scale 
parameter of QCD. At LO, direct photoproduction proceeds via subprocess flT3|) and re¬ 
solved photoproduction via subprocesses (BB and flT^ . The cross section of subprocess 
flT3|) reads [28] 



s 2uQ^ 
f st 



(29) 


dt 


where the Mandelstam variables and gauge coupling constants are dehned in Eqs. (I5|) and 
(]9|), respectively. The ones of subprocesses (ITOj) and (IT^ may be read off from Eqs. (IHHll 
and dSnj), respectively. The emission of photons off the (anti)proton can happen either 
elastically or inelastically, i.e. the (anti)proton stays intact or is destroyed, respectively. In 
both cases, an appropriate PDF can be evaluated in the Weizsacker-Williams approxima¬ 
tion [2^3011311132] . Since these PDFs are of 0{a), the LO photoproduction contributions 
are of 0{a^). Although photoproduction contributions are parametrically suppressed by 
a factor of a/ag relative to the 0 {a'^as) corrections discussed above, we include them in 
our analysis because they may turn out to be sizable in certain regions of phase space. 

We generated the Feynman diagrams using the program package DIANA [33] and 
checked the output using the program package FeynArts 3 [M]. We reduced the one- 
loop tensor integrals to scalar ones using custom-made routines written with the symbolic 
manipulation program FORM version 4.0 [3S]. We evaluated the scalar one-loop integrals 
using the analytic results listed in Ref. [27] . 

3 Numerical analysis 

We are now in a position to present our numerical analysis. As input we use the pole 
masses Mw = 80.385 GeV, Mz = 91.1876 GeV, Mh = 125 GeV, mb = 4.89 GeV, 
and mt = 173.07 GeV, and the MS coupling constants a^M^) = 1/127.944 [36] and 
= 0.1180 [32] to gauge d(/i^) and We employ the NNPDF2.3QED 

NLO set of proton PDFs [32], which also include QED evolution and provide a photon 
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distribution function. This allows us to consistently treat direct and resolved photopro¬ 
duction via inelastic scattering off the (anti)proton along with ordinary hadroproduction. 
In Ref. [32], the QED sector is treated at LO, or, more accurately, at the leading loga¬ 
rithmic level, where factorization is still trivial and does not yet require the specihcation 
of a scheme. In this sense, the NNPDF2.3QED NLO PDFs are compatible with our con¬ 
vention of employing the MS factorization scheme in the QED sector [3^. By the same 
token, the dependence on the QED factorization scheme contributes to the theoretical 
uncertainty, which we refrain from assessing here. As for photoproduction via elastic 
scattering off the (anti)proton, we adopt the photon flux function from Ref. [2S] and the 
resolved-photon PDFs from Ref. |3H|- For dehniteness, we identify the renormalization 
and factorization scales with the Z-boson transverse mass, = \/Pt ^z- 

In Fig. [HI we study the cross section oi pp ^ Z + X at center-of-mass energy \fS = 
1.96 TeV appropriate for Tevatron run II (a) differential in pr integrated over y and 
(b) differential in y imposing the acceptance cut pr > 10 GeV. Specihcally, we show (i) 
the NLO QCD result considered in Ref. [T2|, i.e. the sum of the 0{aas) and O^aa"^) 
results (thin solid lines); (ii) the 0{a^) Born result (thin dot-dashed lines); (iii) the 
purely weak 0{a'^as) corrections to subprocesses fiTOj) and (IT^ considered in Ref. [T7| 
(thick dashed green lines); (iv) the residual electroweak 0 (a‘^as) corrections (thin dashed 
blue lines); (v) the 0 {a^) photoproduction contributions (thin dotted blue lines); and 
(vi) the total sum (thick solid red lines). The pr distributions in Fig. [8](a) are plotted on 
a logarithmic scale. Since the one of contribution (iii) is negative in the considered pr 
range, its modulus is shown. The y distributions in Fig. |8](b) are plotted on a linear scale. 
For better visibility, contributions (ii), (iv), and (v) are amplihed by a factor of 100. In 
Fig. mb), we do not consider negative y values because the y distributions are symmetric 
by charge conjugation invariance. In Fig. [9], we decompose contribution (iv) (thick solid 
red lines) into the combination of the 0{a) QED corrections to subprocesses (ITOD and 
(IT^ and the 0{as) QCD corrections to subprocess (ITTD (thin solid lines), which cannot be 
usefully separated, and the QCD-electroweak interference contributions from subprocesses 
(ITT|) (thin dot-dashed lines) and (IT5]) (thin dashed green lines). In Figs. [TUI and [TTl we 
repeat the analyses of Figs. [H] and [HI respectively, for pp —)■ Z + X at \/S = 14 TeV 
appropriate for the LHC. 

From Figs. [8] and HOl we observe that the combined effect of the electroweak contri¬ 
butions (ii)-(v) is to reduce the NLO QCD predictions (i). The reduction ranges from a 
few percent at low pr values to a few tens of percent in the large-p^ domain, where large 
Sudakov logarithms dominate. The bulk of the electroweak contributions (ii)-(v) is made 
up by the purely weak 0 {a‘^as) corrections to subprocesses fiTOj) and flT^ [contribution 
(iii)], which are negative throughout the kinematic ranges considered here. The residual 
types of electroweak effects taken into accout here, namely the (T(a^) Born result (ii), the 
residual electroweak 0 {a'^as) corrections (iv), and the 0 {a^) photoproduction contribu¬ 
tions (v), are all positive, but numerically suppressed by typically one order of magnitude 
or more relative to contribution (iii), except for contribution (ii) in the small-pT range. 

In Fig. ^Fe QED-type correction is throughout positive, the qq interference con¬ 
tribution is throughout negative, and the qq interference contribution is negative for 
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Pt^ 165 GeV. On the other hand, in Fig. [lira), the QED-type correction is negative 
for GeV, the qq interference contribution is again throughout negative, and the 

qq interference contribution is throughout positive. From Figs, lira) and fTTT a). we observe 
that those two interference contributions are suppressed relative to the QED-type correc¬ 
tion in the small-pr range. This is also reflected in the y distributions of Figs. (HTb) and 
[rTT bi. which receive dominant contributions from the small-pr ranges. Such a suppression 
is expected from the comparison of color factors PBS] . However, the situation is quite 
different at large pt values. In fact, in Fig. |9])a), the qq interference contribution steadily 
approaches the QED-type contribution for increasing value of px, and, in Fig. fTlTa). the 
qq interference contribution exceeds the QED-type contribution for GeV. Gom- 

paring Figs. 191 and ITT] with Figs. IHl and flOl we observe that the three 0{a'^as) contributions 
of class (iv) range at the permille level with respect to the well-known NLO QGD result 
[T2] . Specifically, in the case of a/dy ai y = 0, the QED-type, qq interference, and qq 
interference contributions normalized to the NLO QGD result approximately amount to 

3 X 10“^, 2 X 10““^, and —1 x 10“^ at the Tevatron and to 1 x 10“^, 8 x 10“®, and —2 x 10“^ 
at the LHG, respectively. 

4 Conclusions 

We studied the inclusive hadroproduction of large-pr single Z bosons including both 
the QGD and electroweak NLO corrections and presented px and y distributions under 
Tevatron and LHG kinematic conditions. Our analytic results agree with the literature 
as far as the latter goes. Specifically, we recovered the well-known NLO QGD corrections 
[12], of absolute order and the purely weak one-loop corrections, of absolute order 

0{a‘^as), to the partonic subprocess qq —)■ Zg and its crossed versions [I?]. We completed 
our knowledge of the 0{a‘^as) corrections by providing also the QED corrections and some 
missing weak contributions due to interferences of tree-level 2 —)■ 3 scattering amplitudes 
in compact analytic form ready to be used by the interested reader. While the new 
0 {a'^as) contributions turned out to be numerically small, their knowledge should help 
us to reduce the theoretical uncertainty on this important benchmark cross section. We 
also considered, for the first time, direct and resolved photoproduction in elastic and 
inelastic scattering. 

In the experimental analyses to be compared with our theoretical predictions, Z bosons 
must be identified, preferrably via their decays to e+e“ or pairs, and their four- 

momenta must be reconstructed and sampled in bins of px or y, ignoring any other 
available information about the selected events. If there are more than one identified and 
reconstructed Z boson in an event, then each of them is taken to generate one entry in the 
considered histogram. Theoretical predictions for the hadroproduction of more exclusive 
final states, such as [21] or uu+j [22], require a different mode of experimental data 

analysis. In this sense, the results presented here are not already included in Refs. EM, 
the more so as QGD-electroweak interference contributions of the type mentioned above 
were neglected there. Our detailed study confirmed the expectation EM based on 
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the analysis of + j hadroproduction [20] that those interference contribntions are 
numerically small, for the small-pr regime. 
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A Cross section oiq-\-q^Z-\-X through 0{cx^cXs) 


In this appendix, we consider 0{a‘^as) contributions to q + q^Z + X. This includes the 
virtual QED corrections to subprocess the virtual QGD corrections to subprocess 
dm), the real corrections from subprocess ffT7|) . and the interference contributions from 
subprocess ffTT|) involving a virtual photon or Z boson. We have 
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B Cross section oiq-\-g^Z-\-X through 0{(y.^cXs) 

In this appendix, we consider O(a'^as) contributions to q + g^Z + X. This includes 
the virtual QED corrections to subprocess flT^ and the real corrections from subprocess 
(ED, We have 
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C Cross section oiq-\-q^Z-\-X through 0((x^(Xs) 

In this appendix, we consider 0{a^as) contribntions to q + q^Z + X. This includes the 
interference contributions from subprocess ffTSj) involving a virtual photon or Z boson. 
We have 
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Q2 + S 

Here, we have again used the phase-space integrals given in Eq. (EZD- 


/4(l,2,t, —1) > -|- It-^u 


(45) 
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(a 1) 





Figure 3: Weak-interaction one-loop diagrams contributing to the partonic subprocess 
q + q^Z + gaX 0{a‘^as). The weak-interaction one-loop diagrams contributing to the 
partonic subprocess q + g^Z + q at 0 {a'^as) are obtained by appropriately crossing 
external legs. 
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Figure 4: Tree-level diagrams contributing to the partonic subprocess q^-q^Z + q + q 
at 0{aa‘^) (with a virtual gluon). Interferences of diagrams (a 5)-(a 8) with a virtual 
photon or Z boson (gluon) with diagrams (a l)-(a 4) with a virtual gluon (photon or Z 
boson) contribute at O(a^as). 



Figure 5: Tree-level diagrams contributing to the partonic subprocess q + q^Z + g + g 
at C>(aa^) and to the partonic subprocess q + q^Z + g + j at 0{a‘^as). 
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Figure 6: Tree-level diagrams contributing to the partonic subprocesses q + q ^ Z + q + q 
and q + q' ^ Z + q + q' eit 0{aal). Interferences of diagrams (a 5)-(a 8) with the gluon 
replaced by a photon or Z boson with diagrams (a l)-(a 4) contribute at 0{a'^as). 



Figure 7: 0{a^^'^as) and tree-level diagrams interfering to yield 0{a‘^as) contri¬ 

butions to subprocesses q + q^Z + q + q and q + q—>Z + q + q. 
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(a) 


pp Vs = 1.96 TeV, pt > 10 GeV 



(b) 


Figure 8: Cross section distributions in (a) pt and (b) y for px > 10 GeV oi pp ^ Z + X 
at y/S = 1.96 TeV (Tevatron run II). In each frame, the NLO QCD result [12], i.e. the 
sum of the O(aas) and 0{aa‘^) results (thin solid lines), the (!l(a^) Born result (thin 
dot-dashed lines), the purely weak 0{a‘^as) corrections to subprocesses flTU]) and flT^ [T7] 
(thick dashed green lines), the residual electroweak 0{a‘^as) corrections (thin dashed blue 
lines), the 0{a^) photoproduction contributions (thin dotted blue lines), and the total 
sum (thick solid red lines) are shown. 







(a) 



(b) 


Figure 9: The residual electroweak 0{a‘^as) corrections in Fig. |8] (thick solid red lines) are 
decomposed into the combination of the 0(a) QED corrections to subprocesses (ITOl) and 
(IT^ and the 0{as) QCD corrections to subprocess (ITTl) (thin solid lines), and the QCD- 
electroweak interference contributions from subprocesses (1141) (thin dot-dashed lines) and 
firsj) (thin dashed green lines). 
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(a) 


pp^/S = U TcV, p, > 10 GeV 



(b) 


Figure 10: Same as in Fig. |H1 but ioi pp ^ Z + X at y/S = 14 TeV (LHC). 
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(a) 


pp^/S = U TcV, p, > 10 GeV 



(b) 


Figure 11: Same as in Fig. IHl but ioi pp ^ Z + X at y/S = 14 TeV (LHC). 
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